First-principles study of native point defects in Bi 2 Se 3 
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Using first-principles method within the framework of the density functional theory, we study 
the influence of native point defect on the structural and electronic properties of Bi2Se3. Se 
vacancy in Bi2Se3 is a double donor, and Bi vacancy is a triple acceptor. Se antisite (Sesj) is 
always an active donor in the system because its donor level (e(+l/0)) enters into the conduction 
band. Interestingly, Bi antisite(Bis e i) in Bi2Se3 is an amphoteric dopant, acting as a donor when 
/i e <0.119eV (the material is typical p-type) and as an acceptor when ^i e >0.251eV (the material is 
typical n-type). The formation energies under different growth environments (such as Bi-rich or 
Se-rich) indicate that under Se-rich condition, Sesi is the most stable native defect independent 
of electron chemical potential /i e . Under Bi-rich condition, Se vacancy is the most stable native 
defect except for under the growth window as /i e >0.262eV (the material is typical n-type) and 
A/xse<-0.459eV(Bi-rich), under such growth windows one negative charged Bis e i is the most stable 
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I. INTRODUCTION 



The narrow-band-gap semiconductor Bi2Se3 
(E g ~0.35eV)i2 has been best known for a long 
time as an excellent thermoelectric material because 
of their unique near-gap electronic structure and high 
thermoelectric figure of merit ZT=S 2 aT/Ki2r— Recently, 
with the development of topological insulators (TIs),— ~— 
it attracts the research attention again. TI is a new state 
having an energy gap in its bulk band structure and 
metallic helical states on its surface, which distinct from 
simple metal or insulator Bi2Se3 is a strong three 
dimensional (3D) topological insulator with its surface 
states consisting of single Dirac cone at the T point 
which is protected by the time-reversal symmetry from 
any time-reversal perturbation, such as crystal defects 
and nonmagnetic impurities^ ' 11 ' 12 Particularly, its bulk 
band gap is much larger than the energy scale of room 
temperature, making it the most suitable candidate for 
the high-temperature spintronics applicationi&ii 

Whether as thermoelectric material or topological 
insulators, the native defects in Bi2Se3 play important 
role in the material. As thermoelectric material, it needs 
well-defined electrical and thermal conductivities, high 
mobility of free current carriers, and thermoelectric 
power, all of them will be influenced by the presence 
of native defectsJ^r— As a topological insulator, the 
intrinsic defects such as antisites or vacancies behave 
as n type dopant and consequently shift the Fermi 
level above the Dirac point, which makes it difficult to 
characterize the topological transport properties and to 
realize topological devices because both of them severely 
rely on the behavior of surface Dirac fermions i 12 ' 16 ' 17 
Hence, it is vital to investigate the nature of native 
defects in Bi2Se3. In experiments, lots of effort o 13 ' 18 



have been pursued on such issue. Urazhdin et alj^ 3 - 
have found that Bi2Se3 single-crystal behaves as n 
type material under excess Bi or Se growth condition. 
They predicted that Bis e antisites contribute shallow 
acceptors when the samples are doped with excess Bi. 
However, the solubility of Bi in Bi2Se3 is very low and 
the n type doping behavior of the samples is derived 
from the compensating defects. While when the samples 
doped with excess Se, Sesi antisites introduce shallow 
donor in the materials. 13 Bludska et alJ^ concluded 
that the presence of overstoichiometric Bi atoms can 
induce the following defects: (i) Se vacancy carrying two 
positive charges (Vjt 2 ), which is attributed to the main 
reason for n-type Bi2Se3; (ii) Bis e antisite carrying one 
negative charge (Bi^ 1 ); (iii) Bi vacancy supposed to 
carry three negative charges (V^ 3 ); (iv) other complex 
defects. In theoretical calculations, Shuang-Xi Wang 
et alJ^ have studied the formation energies of native 
point defects in Bi2Se3, but the study only gave the 
results under two extreme growing conditions: Se-rich 
and Bi-rich. The systematic study of the evolution of 
the formation energies according to growth conditions 
has not been reported as yet. 

To this end, in our present work, the stability and 
formation energies of native defects in Bi2Se3 are 
systematically investigated using first-principles method 
based on density function theory. The results indicate 
that Sel vacancy is a double donor, Bi vacancy is a 
triple acceptor, and Sesi is an always active n-type 
dopant in Bi2Se3. Interestingly, we find that Bis e i is 
an amphoteric dopant. By analyzing the formation 
energies of all defects under different growth conditions 
(such as Bi-rich or Se-rich), we find that under Se-rich 
growth condition, Se^i is the most stable one which is 
independent of electron chemical potential [i e . However, 
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FIG. 2: Locations of three types of vacancies (a) and two 
types of antisite (b) defects. 



FIG. 1: Rhombohedral primitive (a) and the hexagonal con- 
ventional unit cells (b) of Bi2Se3. 



TABLE I: Lattice parameters-^— in our calculation for 
Bi(space group R3m), Se(space group P3i21), and Bi2Se3. 



System 


Ref. 


a(A) 


c/a 






Bi 2 Se 3 


28 


4.138 


6.92 


0.399 


0.206 


Bi 


29 


4.535 


2.611 


0.234 




Se 


30 


4.368 


1.135 




0.225 



under Bi-rich growth condition, Se vacancy is the most 
stable defect, except for a small growth window which is 
defined by /i e >0.262eV (the material is typical n-type) 
and A/us e <-0.459eV(Bi-rich). In that growth window, 
one negative charged Bis e i is the most stable native 
defect. 



II. COMPUTATIONAL METHOD AND 
PROCEDURES 

Bi2Se3 has the rhombohedral crystal structure with 
space group D| d (R3m). Its primitive unit cell contains 
three Se atoms and two Bi atoms, and Se atoms have two 
nonequivalent sites Sel and Se2, as shown in Fig. 1(a). 
The conventional hexagonal unit cell has a layered struc- 
ture with three basic unit cells (one unit cell named a 
quintuple layer(QL)) weakly bound to each other by the 
van der Waals forces i 20 ' 21 In each QL, five atomic planes 
with atomic order Sel-Bi-Se2-Bi-Sel, see as fig. 1(b), and 
the chemical bond between Bi and Se atoms is of the 
covalent-ionic type^ Our calculations are performed in a 
3x3x1 supercell with 135 atoms by using the plane- wave 
pseudopotential code Vienna ab initio simulation package 



(VASP)i 23 ' 24 The core-electron interaction is modeled by 
the projector augmented wave metho d 25 ' 26 and the gen- 
eralized gradient approximation of Perdew, Burke, and 
Ernzerhof, exchange-correlation functional 2 ^ is adopted. 
The plane wave cutoff energy is chosen as 240eV, and 
the Brillouin zone is sampled by using 5x5x1 Gamma- 
centered Mon-khorst-Pack grids. All structure parame- 
ters of Bi2Se3 are chosen from experimental data as listed 
in Table I. The energy convergent criterion is 10 eV per 
unit cell. We relax all atoms without spin-orbit coupling 
(SOC) because the SOC only slightly affect the relax- 
ation results and is too much time consuming. The forces 
on all relaxed atoms are less than O.OleV/A. The self- 
consistent calculations are performed with and without 
spin-orbit coupling (SOC) to make a comparation. We 
find that the formation energies with SOC are lower than 
that without SOC and the difference of transition energy 
levels is larger than 0.04eV. In this paper, we emphasize 
on the results with SOC. Our calculation gives a band 
gap of 0.31eV, which is in good agreement with experi- 
mental and other theoretical reportsi 2 ' 11 ' 31 ' 32 

We consider three types of vacancy point defect in 
Bi2Se3, namely vacancy on the Sel, Se2 and Bi sublat- 
tices, denoted as Vs e ii Vse2, and N Bi- We sign their 
detail locations in Fig. 2 (a). As for antisite defects, we 
consider Se^ (substitute one Bi atom by Se atom) and 
Bis e i (substitute one Sel atom by Bi atom) because ac- 
cording to the formation energies discussed below, the 
formation energy of Vse2 is 0.4eV higher than that of 
Vs e i, as shown in Fig. 2 (b). 

It is well known that the calculations based on GGA 
typically underestimate the band gap of a semiconductor 
and cannot give the absolute position of the defect level. 
In this paper, we emphasize the analysis of the defect for- 
mation energy and the transition energy levels! 33 ' The 
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FIG. 3: (a)Formation energies of all defects as a function of the electron chemical potential. The slope of the lines gives 
the charge state. (b)Formation energy of each defect under different growth condition, such as Se-rich(A/is e =OeV) and Bi- 
rich(A/isi=OeV or A^ts e =-0.491eV). (c) is the project view figure of figure (b). The VBM of all systems are set to zero. 



formation energies of all defects in Bi2Se3 are given by: 

rr(defect,q) _ rj(defect,q) ^perfect a 
H f - ^totle ~ ^totle ~2^ An i^ 

i 

+ 1 E VBM + 9Me (1) 

where }{ye-fect,q) ^ e formation energy of a defect in 
Bi 2 Se3 under charge state q. The terms, E^/ e ect ' q ' and 

Etotil eCt ' are the total energy of the super-cell with de- 
fect under charge state q and pure Bi2Se3 in the same 
super-cell size, respectively. An^ is the change in the 
number of atoms of species i, which represents Se or Bi 
atom in the present study, /ii is the chemical potential of 
element i. E v e g^f is the valence band maximum (VBM) 
of the super-cell with defect. /x e is the electron chemical 
potential. The defect transition energy level e(q/q') is 
the electron chemical potential in Eq. (1), at which the 
formation energy of the defect with charge q is equal to 
that of another charge q'. 

To describe formation energy more accurately, we in- 
troduce two corrections as follows: i) The total energies 
of the charged systems should be corrected for the inter- 
action of the charged defect with the compensating back- 
ground and its periodic images. We use Makov-Payne 
(M-P) corrections, 35 formulated as q 2 a/2eL, where L is 
the linear dimension of the supercell, e is the static di- 
electric constant and a is Madelung constant; ii) Because 
the defects induce significant distortion of the band struc- 
ture and fluctuation of the bandgap, the Evbm obtained 
from the first-principles method directly can not be used 
in Eq.(l). We adopt the correction introduced by previ- 
ous works4^— Firstly, we assume that the potentials in 
the perfect super-cell are similar to those far from a defect 
in a defective super-cell. Then, the average potential of 
the plane farthest from the defect in the defective system 
{Yav ) anc l the average potential of the corresponding 
plane in the perfect system (V^ r ^ ect ) are determined. 
The difference of the average potentials between the per- 



fect and defective super-cells is used to determine Evbm 
of the defective super-cell as follows: 

rpdefect _ ^perfect , ^defect _ rrperfect 
VBM ~ VBM t v av v au l z i 

The first term on the right-hand side of Eq. (2) is 
the VBM of perfect super-cell and can be obtained by 
E v % f ^ ct = E T (per feet : 0) - E T (per feet : +1), where 
Et (perfect : q) is the total energy of a perfect super-cell 
under charge state q. 



III. RESULTS AND DISCUSSIONS 
A. The transition energy levels 

We investigate the transition energy levels of all 
defects quantitatively by calculating their formation 
energies in function of the electron chemical potential 
/i e , and the chemical potentials /x; are considered to 
be the total energies of ideal elemental bulk Bi and Se 
crystal simply(such choice does not influence the results 
of the transition energy level of native defect), as shown 
in Fig. 3 (a). For Se vacancy, we consider both Vs e i and 
Vse2, and find that the formation energy of Vg e 2 is 
0.4eV higher than that of Vs e i- The lower formation 
energy of Vsei can be explained as a result of weak 
interaction between QLs. Our results indicate that the 
stable charge states of Vs e i in Bi 2 Se 3 are Vg^, Vg el 
and Vg el as the Fermi level moves through the bandgap, 
which shows that Vs e i acts as a double donor in the 
systems. The results is in good agreement with the 
previous report^ The donor levels e(+2/+l)=0.209eV, 
e(+l/0)=0.231eV, are O.lOleV and 0.079eV below the 
conduction band minimum(CBM) respectively. As for 
Sesi in Bi2Se3, it is one positive charged as /i e moves 
through the bandgap energy window. Consequently, 
Sesi in Bi2Se3 behaves as a donor, and the donor level 
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e(+l/0) enter into the conduction band, indicating that 
it is an active n-type doping in the system. The anti-site 
Se protrudes outwards about 0.14A to the Sel layer and 
induce the distance between the Se atoms in the Sel 
layer decrease about 10.56%. The stable charge states 
of isolated Bi vacancy (Vsi) in Bi2Se3 as the Fermi 
level moves through the bandgap are VlL, V B ^ and 
V^f , as shown in Fig. 3(a), indicating its triple acceptors 
characteristics. The acceptor levels are e(0/-l)=-0.07eV 
(the level is active in the system because it is below 
the valence band maximum (VBM) ), e(-l/-2)=0.008eV 
(shallow acceptor level), and e(-2/-3)=0.129eV (deep 
acceptor level), respectively. The results agree with the 
previous report. 18 The last native defect considered in 
our present work is Bisei- The results of formation 
energy of Bis e i as shown in Fig. 3(a) indicates that 
the stable charge states of isolated Bis e i in Bi 2 Se 3 
as the Fermi level moves through the bandgap are 
-^Sei' Bi^ei and Big el . Its transition energy levels are 
e (+l/0)=0.119eV and e(0/-l)=0.251eV. This means 
that the Bis e i in Bi2Se3 is an amphoteric dopant, acting 
as a deep donor when the electron chemical potential 
is close to VBM (the material is typical p-type) and as 
a deep acceptor when the electron chemical potential 
is close to CBM (the material is typical n-type). Bis e i 
will behave as a compensation dopant to the electrical 
type of the material. The result agrees well with the 
band structure of the Bi 2 Se3(lll) surface with defect 
Bis e i previously reported^ with no apparent donor or 
acceptor feature. Moreover, in experiment, the n-type 
Bisei was not mentioned and we will discuss the reasons 
below. As for the geometric configuration around Bis e i, 
we find that the anti-site Bi atom protrudes outward 
about 0.419A, 0.425A and 0.439A relative to the average 
height of Sel layer under neutral, one negative and one 
positive charge state, respectively. Such results are in 
good agreement with the experimental observation^, 
that there are protrusion point on the surface when 
Bi2Se3 doped with excess Bi. The protrusion can be 
attributed to the anti-site Bi on the surface denoted as 
Sei layer in the present work. 



TABLE II: Formation energies (eV) of each defect at typical 
Se-rich (A/iSe=0eV) and typical Bi-rich(A/iSe=-0.491eV or 
AfiBi~0) growth conditions while ^i e =0eV and /^ e =0.31eV. 



Afi Se =0eV 
System /j, e —0eV /i e =0.31eV 
Vsei 0.704 1.144 
V Bi 1.867 1.074 

Sesi 0.529 0.839 
Bisei 1.773 1.833 



A/i S i=0eV 
^ e =0eV /i e =0.31eV 
0.213 0.653 
2.604 1.811 
1.756 2.066 
0.545 0.605 



Bi2Se3 satisfy the relationship: 



(3) 



where Msi^Sea * s the chemical potential of one formula 
unit of Bi2Se3. The chemical potentials of Bi and Se in 
ideal elemental bulk are denoted as n° Bi and Hg e , respec- 
tively. Then the the difference of the chemical potential 
ofBi(Se) is defined as: Afi Bi =fi Bi -fi Bi (A^ Se =/i Se -A<s e ). 
The Gibbs free energy in the process of formation of one 
formula unit of Bi2Se3 is given as: 



AG(Bi 2 Se 3 ) = n» B tse 3 
= 2An Bi 



- 2»° m - 3/4 e 

- 3A/i Se 



(4) 



From Eq. [TJ [3] and 2J the formation energy formula 
of point defects in Bi2Se3 as a function of Afise can be 
described as: 
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(def ect,q) 



-AGAriBi - {Arise - - An Bi )AMs e (5) 
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(defect,q) 
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^perfect 
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Similarly, the formation energy formula as a function of 
A/i B i also can be obtained as above. Eq.[5]and[6]indicate 
that the formation energy is a function of [i e and A[ig e . 
In addition, Eq.(4) tells us that to form Bi 2 Se3, Afi B i 
and AjjLSe should meet the following condition: 



B. Stability of the native defect 



1 



-AG < Afi Bl < 



(7) 



We analyze the stability of all the native defects 
by calculating their formation energies under different 
growth environments (such as Bi-rich or Se-rich) us- 
ing the method proposed by Hashibon et al4S Using 
Afj,Bi(Afj,Se), the difference of the chemical potential of 
Bi(Se) between in the bulk Bi2Se3 and in ideal elemental 
bulk Bi(Se), the method can accurately simulate the re- 
lationship between growth condition and the formation 
energy of native defect. The detail process of this method 
is as follows: 

Firstly, the chemical potential of Bi and Se in bulk 



-AG < An Se < (8) 

The formation energies under different grown envi- 
ronments are shown in Fig. EJb) and the typical results 
for /i e =0eV, ^ e =0.31eV under A/xs e =0eV (Se-rich) and 
AfiBi—0eV (Bi-rich) conditions are listed in Tab. [11] The 
results as shown in Fig. [3jb) and (c) indicate that under 
Se-rich growth condition, the growth window denote as 
I in Fig.|3Kc), SeBi is the most favorable native defect in 
the system within the full energy window of the electron 
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FIG. 4: Total DOS of perfect Bi,g e system and Bis e system with native point Vsei, SeBi and Bis e i respectively. The Fermi 
level is set to zero and denoted by the vertical dashed line. 



chemical potential /J, e , it behaves as a single donor. This 
means that under Se-rich growth condition Sesi is the 
predominant native defect in the system and responsible 
for the n-type characteristics of the system under such 
growth condition. For example, when A/Lts e =0eV and 
/i e =0 eV, the formation energy of Se^i is 0.175eV, 
1.338eV and 1.244eV lower than that of V Se i, V Sl , and 
Bis e i, respectively. When /^s e =0eV and /^ e =0.31eV, 
the formation energy of Sent is 0.305eV, 0.235eV, and 
0.994eV lower than that of Vsei, V_Bi, and Bis e i, 
respectively. When the growth window is in the region 
II as shown in Fig.(3jc), the most favorable native defect 
is Vs e i, whereas in the growth window III which is 
confined by // e >0.262eV (the material is typical n-type) 
and A/j,s e <-0.459eV(Bi-rich) one negative charged 
Bisei is the most favorable one. This means that under 
Bi-rich growth condition, Vs e i is the predominant native 
defect and it is the main reason for n-type electrical 
characteristics of the systems. For example, when 
AfiBi—OeV (A/is e =-0.491eV), the formation energy of 
Vs e i according to /i e =0 eV is 0.213eV, which is 2.391eV, 
1.543eV and 0.333eV lower than that of Vsi, Sesi and 
Bis e i, respectively. Only when A^se<-0.459eV (the 
material is typical n-type), p-type Big e i is the dominant 
native defect and it behaves as a compensation dopant 
for the the material. The formation energy of Big e i 
according to /i e =0.31eV is 0.605eV, which is 0.048eV, 
1.206eV and 1.461eV lower than that of Vsei, V b% and 
Sesi, respectively. 

In experiments, the electronic properties of bulk 
Bi2Se3 crystals are usually dominated by electron 
donors, resulting in n-type conductivity. Based on 
above analysis, such n-type conductivity mainly derives 
from two reasons: First, the formation energy of p-type 



native defect Vbi is higher than that of n-type native 
defect, such as Se^i and Vs e i. Second, one negative 
charged p-type native defect Bigei is only stable when 
the material is strong n-type and Bi-rich growth con- 
dition (A/X5 e <-0.459eV) behaving as a compensation. 
Moreover, we also find that one positive charged Bigei 
only forms when material is p-type and its formation 
energy is higher than that of Vsei which is difficult to 
be found in experiment. 



IV. ELECTRONIC STRUCTURES 

The density of states (DOS) of perfect Bi 2 Se3 system 
and three systems with native defect of Vsei, Se^i or 
Bisei are shown in Fig. 2] In comparison with that of 
perfect system, it is clearly seen that each defect induces 
two defect states into the band gap of the material. For 
Vsei, the defect states just below the Fermi level mainly 
come from p states of three nearest neighbor (NN) Bi 
atoms of Vs e i and are completely filled by two electrons. 
The electrons filled the defect states are easily released 
leading the defect a double donor. The defect states of 
Sesi mainly come from p states of four Se atoms (the Se 
atom substituted the Bi and its three NN Se atoms), 
which arc half filled. The electron filled the defect state 
is easily released leading the defect a single donor. The 
defect states of Bis e i mainly come from p states of four 
Bi atoms (the Bi atom substituted Sel and its three 
NN Bi atoms), which are half filled. Depending on the 
conductive type of the material, either the electron filled 
the defect states is released leading the defect a donor or 



6 



the unfilled state capture one electron leading the defect 
an acceptor. Bis el behaves as an amphoteric dopant in 
the system. 

Comparing the DOS and the formation energy results, 
we find that the defect levels in DOS are lower than the 
results of formation energy. The phenomenon may come 
from two possible reasons as follows: (i) the different 
relaxed geometries between neutral and charged system. 
Taking Vs e for example, the distance between Bi 
atoms (nearest to the Sel layer with vacancy) changes 
obviously from 3.834A for q=0 to 4.370A for q=+2; 
(ii)the correction of VBM for each charged system tend 
to make both acceptor and donor level higher. Moreover, 
we also note that although both Sesi and Bisd induce 
half filled defect states around the Fermi level, Sesi acts 
as a single donor whereas Big e i behaves as amphoteric 
dopant. The reason is that the un-filled defect level 
above the Fermi level of Sest is higher than that of the 
Bis e i, which make it difficult to capture electron acting 
as p-type dopant. 



V. CONCLUSION 

In present work, we find that Vs e in Bi2Se3 is a double 
donor; V Bi is a triple acceptor; Sesi is an active donor, 



whereas Bigei is an amphoteric dopant acting as a donor 
when /Li e <0.119eV(the material is typical p-type) and as 
an acceptor when /j e >0.251eV(the material is typical 
n-type) . The analysis of the stabilities of all defects show 
that under Se-rich growth condition, Se^ is most stable 
according to all electron chemical potential fi e . Under 
Bi-rich growth condition, Vs e is the main defect. One 
negative charged Risei is most stable in a small growth 
window which defined by /x e >0.262eV (the material is 
typical n-type) and A^5 e <-0.459eV(Bi-rich). 
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